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Abstract 

The microcanonical multifragmentation model from [Al. H. Raduta and Ad. 
R. Raduta, Phys. Rev. C 55, 1344 (1997); 56, 2059 (1997); 59, 323 (1999)] is 
refined and improved by taking into account the experimental discrete levels 
for fragments with yl < 6 and by including the stage of sequential decay of 
the primary excited fragments. The caloric curve is reevaluated and the heat 
capacity at constant volume curve is represented as a function of excitation 
energy and temperature. The sequence of equilibrated sources formed in the 
reactions studied by the ALADIN group (^*^'^Au+^*^'^Au at 600, 800 and 1000 
MeV/nucleon bombarding energy) is deduced by fitting simultaneously the 
model predicted mean multiplicity of intermediate mass fragments (Mjmf) 
and charge asymmetry of the two largest fragments (012) versus bound charge 
inbound) on the corresponding experimental data. Calculated HeLi isotopic 
temperature curves as a function of the bound charge are compared with the 
experimentally deduced ones. 
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I. INTRODUCTION 



Nuclear multifragmentation is presently intensely studied both theoretically and exper- 
imentally. Due to the similitude existent between the nucleon-nucleon interaction and the 
van der Waals forces, signs of a liquid-gas phase transition in nuclear matter are searched. 
While the theoretical calculations concerning this problem started at the beginning of 1980 

the first experimental evaluation of the nuclear caloric curve was reported in 1995 by 
the ALADIN group 0. A wide plateau situated at around 5 MeV temperature lasting from 
3 to 10 MeV/nucleon excitation energy was identified. The fact was obviously associated 
with the possible existence of a liquid-gas phase transition in nuclear matter and generated 
new motivations for further theoretical and experimental work. Similar experiments of EOS 
0] and INDRA followed shortly. Using different reactions they obtained slightly differ- 
ent caloric curves, the plateau - like region being absent in the majority of cases. Factors 
contributing to these discrepancies are both the precision of the experimental measurements 
and the finite-size effects of the caloric curve manifested through the dependency of the 
equilibrated sources [-E'*(A)] sequence on the reaction type. 

Concerning the first point of view, recent reevaluations of the ALADIN group concerning 
the kinetic energies of the emitted neutrons brought corrections of about 10 % (in the case of 
the reaction ^^'^Au+^^''Au, 600 MeV/nucleon). More importantly however it was proven that 
the energies of the spectator parts are growing with approximately 30 % in the bombarding 
energy interval 600 to 1000 MeV/nucleon. On the other side, the universality of the quantity 
MiMpiZbound) subject to the bombarding energy variation (which was theoretically proven 
to be a signature of statistical equilibrium) suggests that for the above-mentioned 
reactions the equilibrated sources sequence [i?*(A)] should be the same. Consequently, we 
deal with an important nonequilibrium part included in the measured source excitation 
energies which may belong to both pre-equilibrium or pre-break-up stages [0. The SMM 
calculations suggest a significant quantity of nonequilibrium energy even in the case of the 
600 MeV/nucleon bombarding energy reaction 
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Thus, the necessity of accurate theoretical descriptions of the break-up stage and of the 
sequential secondary particle emission appears to be imperative in order to distinguish be- 
tween the equilibrium and nonequilibrium parts of the measured excitation energies. These 
approaches should strictly obey the constrains of the physical system which, in the case of 
nuclear multifragmentation, are purely microcanonic. As we previously underlined pi]| , P 



in spite of their success in reproducing some experimental data, the two widely used statis- 
tical multifragmentation models (SMM and MMMC |jl3|) are not strictly satisfying the 
microcanonical rules. 

The present paper describes some refinements and improvements brought to the sharp 



microcanonical multifragmentation model proposed in |[T^JT5[] and also the employment of 
the model in its new version in the interpretation of the recent experimental data of the 
ALADIN group MM. 



The improvements brought to the model 1 14,15] are presented in Section II. Section III 



presents the new evaluations of temperature curves and the first evaluations (performed 
with this model) of heat capacities at constant volume {Cy) represented as a function of 
system excitation energy and temperature and also the comparison between the model pre- 
dictions and the recent experimental HeLi isotopic temperature curve [T//eLi(^bo«nd)] 
Conclusions are drawn in Section IV. 

II. IMPROVEMENTS BROUGHT TO THE MICROCANONICAL 
MULTIFRAGMENTATION MODEL 

The improvements brought to the microcanonical multifragmentation model concerns 
both the break-up stage and the secondary particle emission stage, 
(i) Primary break-up refinements 

Comparing to the version of Ref. |l4| , p!5| , p^ the present model has the following new features: 



(a) The experimental discrete energy levels are replacing the level density for fragments with 
A < 6 (in the previous version of the model a Thomas Fermi type level density formula was 
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used for all particle excited states). In this respect, in the statistical weight of a configuration 
and the correction factor formulas [l^Jl^ the level density functions are replaced by the 
degeneracies of the discrete levels, (2Si + 1) (here Si denotes the spin of the ith. excited 
level). As a criterion for level selection (i.e. the level life-time must be greater than the 
typical time of a fragmentation event) we used F < 1 MeV, where F is the width of the 
energy level. 

(b) In the case of the fragments with A> Q the level density formula is modified as to take 
into account the strong decrease of the fragments excited states life-time (reported to the 
standard duration of a fragmentation event) with the increase of their excitation energy. To 
this aim the Thomas Fermi type formula [|1^ is completed with the factor exp(— e/r) (see 
Ref. m): 

P^"^ = ^ exp (2v^) exp (-e/r), (2.1) 
ev48 

where a = A/a, a = 4.7(1.625 + e/B{A, Z)) and r = 9. 
(ii) Inclusion of the secondary decay stage 

For the A > Q nuclei it was observed that the fragments excitation energies are sufficiently 
small such as the sequential evaporation scheme is perfectly applicable. According to Weis- 
skopf theory [|l^] (extended as to account for particles larger than a), the probability of 
emitting a particle j from an excited nucleus is proportional to the quantity: 

' tJo ^-f^' P(E') ^ 

where ej are the stable excited states of the fragment j subject to particle emission (their 
upper limit is generally around 7-8 MeV), E is the kinetic energy of the formed pair 
in the center of mass (cm.) frame, gj = 2Si + 1 is the degeneracy of the level i, fij 
and Bj are respectively the reduced mass of the pair and the separation energy of the 
particle j and finally aj is the inverse reaction cross-section. Due to the specificity of the 
multifragmentation calculations we considered the range of the emitted fragments j up to 
the A = 16 limit. For the inverse reaction cross-section we have used the optical model 



based parametrization from Ref. |T8[. The sequential evaporation process is simulated by 
means of standard Monte Carlo (see for example [0). 

For nuclei with A < A < 6 (the only excited states of A = 4 nuclei taken into consideration 
are few states higher than 20 MeV belonging to the a particle) depending on their amount 
of excitation we consider secondary break-up for e > B{A, Z)/3 and Weisskopf evaporation 
otherwise (here e is the excitation energy of the fragment {A, Z) and B{A, Z) is its binding 



energy). The microcanonical weight formulas have the usual form excepting the level 
density functions which are here replaced by the discrete levels degeneracies. Due to the 
reduced dimensions of the A < 6 systems, the break-up channels are countable (and a 
classical Monte Carlo simulation is appropriate) when a mean field approach is used for the 
Coulomb interaction energy. In this respect, the Wigner-Seitz approach [jl2[ is employed for 
the Coulomb interaction: 

where A and Z denotes the mass and the charge of the source nucleus, the resulting fragments 
have the index i, R%zJ R^z, = (1 + [{Zi/ A,) / {Z / A)f'^ and V = [l + k)Vq. Here V 
denotes the break-up volume and Vq the volume of the nucleus at normal density. It should 
be added that R is the radius of the source nucleus at break-up and -Ra^z, is the radius of 
fragment i at normal density. 

For each event of the primary break-up simulation, the entire chain of evaporation and 
secondary break-up events is Monte Carlo simulated. 



III. RESULTS 

Using the improved version of the microcanonical multifragmentation model, the caloric 
curves corresponding to two freeze-out radii (R=2.25 A^^^ and R=2.50 A^/^ fm) are reevalu- 
ated for the case of the source nucleus (70, 32) (the microcanonical caloric curves evaluated 
with the initial version of the model are given in Ref. pil|JTT| , |2(][| ) . These are presented in 
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Fig. 1 (a). One can observe that the main features of the caloric curve from Refs. pO| , [TT|j20 
are reobtained. Thus, one can recognize the hquid-hke region at the beginning of the caloric 
curve, then a large plateau-like region and finally the linearly increasing gas-like region. 
One may also notice that the caloric curve behavior at the freeze-out radius variation is 
maintained: The decrease of the freeze-out radius leads to a global lifting of the caloric 
curve. 

As it is well known, the curves of the constant volume heat capacity (Cy) as a function of 
system excitation energy (E*) and as a function of temperature (T) may provide important 
information concerning the transition region and the transition order. For this reason the 
curves Cv{E*) and Cv(T) have been evaluated (see Fig. 1 (a) and Fig. 1 (b)). We remind 
that the constant volume heat capacity (Cy) is calculable in the present model using the 



formula 1 1 



Cy' = l- 



3.. 5\ 
-Nr - - — 
2 ^ 2j K_ 



It can be observed that the Cv{E*) curve has a sharp maximum around 4.5 MeV/nucleon 
excitation energy for both considered freeze-out radii. This suggests that a phase transition 
exists in that region. The transition temperatures can be very well distinguished by analyzing 
the Cv{T). One can observe two sharp-peaked maxima pointing the transition temperatures 
corresponding to the two considered freeze-out radii. 

In order to make a direct comparison between the calculated HeLi isotopic temperature 
and the recent experimental results one has to deduce the sequence of excitation energy 
as a function of the system dimension This is done as in Refs. using as match- 

ing criterion the simultaneously reproduction of the (Mjmf) {{Zbound)) and (012) {{Zbound)) 
curves. This couple of curves can fairly well identify the dimension and the excitation of 
the equilibrated nuclear source Here Mjmf stands for the multiplicity of intermediate 
mass fragments and is defined as the number of fragments with 3 < Z < 30 from a fragmen- 
tation event while denotes the charge asymmetry of the two largest fragments and, for 
one fragmentation event is defined as = [Zmax — Z2)/{Zmax + Z2) with Zmax < ^2 < 2 
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where Zmax is the maximum charge of a fragment and Z2 is the second largest charge of a 
fragment in the respective event. Zhound represents the bound charge in one fragmentation 
event and is defined as the sum of the charges of all fragments with Z > 2. 

The simultaneous fit of the calculated curves {Mjmf) {{Zbound)) and (012) ({Zbound)) on 
the corresponding experimental data (^^^Au+^^''Au at 1000 MeV/nucleon) is given in Fig. 
2. The agreement is very good. The equilibrated source sequence we used for 

this purpose is given in Fig. 3 together with the experimental evaluations of the excitation 
energies as a function of source dimension for the reaction ^^''Au+^^^Au at 600, 800 and 1000 
MeV/nucleon. The theoretically obtained sequence is relatively close to the experimental 
line corresponding to 600 MeV/nucleon bombarding energy. The deviations between the 
calculated equilibrated source sequence and the three experimental lines suggest that the 
experimental evaluations contain a quantity of non-equilibrium energy which grows with 
increasing the bombarding energy. As suggested in Ref. |0J^, its origin may be situated in 
both the pre-equilibrium and pre-break-up stage. These deviations are exclusively due to 
the neutron kinetic energies which, reevaluated from the 1995 data [0], are much larger. 

It should also be pointed that apart from the SMM predictions the quantity of non- 
equilibrium energy predicted by the present model is smaller and thus the model predicted 
equilibrated source sequence is closer to the experimental line of the 600 MeV/nucleon 
bombarding energy reaction. 

After evaluating the sequence of the equilibrated sources a direct comparison the HeLi 
calculated isotopic temperature curve with the ones recently evaluated by the ALADIN 
group [^,|[ is performed. To this purpose the uncorrected Albergo temperature is used: 
TneLi = 13.33/ In [2.18 {YeLi/YTu) / (YiHe/YHe)], the experimental predictions being di- 
vided by fx = 1.2 (which is the factor used in the ALADIN evaluation of the HeLi caloric 
curve chosen as to average the QSM, GEMINI and MMMC models predictions). The result is 
represented in Fig. 4 as a function of Zbound- It can be observed that the agreement between 
the calculated TueLi^Z bound) and the experimental data corresponding to the ^^^Au+^^''Au 
reaction at 600 and 1000 MeV/nucleon bombarding energy is excellent on the entire range 
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of Zbound- In comparison, the SMM model predicts in the region Zbound < 25 a curve steeper 
than the experimental data. 

IV. CONCLUSIONS 

Sumarizing, the microcanonical multifragmentation model from Ref. |]n|,|10| is improved 
by refining the primary break-up part and by including the secondary particle emission. 
The caloric curve rededuced with the new version of the model preserves its general aspect 
10| , pT|pO| manifesting an important plateau-like region. The transition regions are clearly 



indicated by the sharp maxima of the Cv{E*) and Cy(T) curves. The model proves the 
ability of simultaneously fitting the "definitory" characteristics of the nuclear multifrag- 
mentation phenomenon (Mjmf) {{Zbound)) and (012) {{Zbound))- Evaluating the equihbrated 
source sequence E*{A) [by using the criterion of reproducing both (Mjmf) and (012) versus 
{Zbound)], a nonequilibrium part of the experimentally evaluated excitation energy growing 
with the increase of the bombarding energy is identified. The direct comparison of the cal- 
culated HeLi caloric curve shows an excellent agreement with the experimental HeLi curves 
recently evaluated by the ALADIN group. 
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FIG. 1. Microcanonical temperature as a function of source excitation energy (a), heat ca- 
pacity at constant volume as a function of source excitation energy (b) heat capacity at constant 
volume as a function of microcanonical temperature (c) . Calculations have been performed for the 
source nucleus (70, 32) with two values of the freeze-out radius. 
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FIG. 2. {Mjmf) {{Zbound)) and (012) {{Zbound)) evaluated by means of the microcanonical 
model in comparison with the experimental evaluations corresponding to the reaction ^^^Au+''^^^Au 
at 1000 MeV/nucleon bombarding energy. (The deviation of calculated (012) from experimental 
data for Z})ound > 65 is, as explained in [7], due to some detection problems.) 
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FIG. 3. The sequence of equilibrated sources evaluated by means of the microcanonical model 
(open circles) following the criterion of simultaneous fitting of the calculated MjMF{Zbound) and 
auiZbound) on the corresponding experimental data. The close symbols represent the sequence 
of excitation energy experimentally measured for the reaction ^^^Au+^^^Au at 600, 700, 1000 
MeV/nucleon bombarding energy [7]. 
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FIG. 4. HeLi temperature curves evaluated with the microcanonical model (open circles) 
in comparison with the experimental HeLi temperatures [7,8] corresponding to the reactions 
^^^Au+^^^Au at 600 and 1000 MeV/nucleon bombarding energy. 
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